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Summary: By means of the linear dynamic conductivity, inductively 
measured on epitaxial films between 30mHz and 30 MHz, the transition 
line Tg{B) to generic superconductivity is studied in fields between B=0 
and 19T. It follows closely the melting line Tm(B) described recently in 
terms of a blowout of thermal vortex loops in clean materials. The critical 
exponents of the correlation length and time near Tg{B), however, seem 
to be dominated by some intrinsic disorder. Columnar defects produced 
by heavy-ion irradiation up to field-equivalent-doses of Btf, — lOT lead to 
a disappointing reduction of Tg {B 0) while for B > B^ the generic line 
of the pristine film is recovered. These novel results are also discussed in 
terms of a loop-driven destruction of generic superconductivity. 



1 Introduction 

Due to the complex electronic and real lattice structure of the high- Tc cuprates 
the nucleation of long-range superconductivity is still heavily debated. In the 
first place, this applies to the pairing mechanism in a non-Fermi liquid in a 
doped Mott-insulator Q giving rise to the formation of a local superfluid den- 
sity near some mean field temperature Tq. Second, the appearance of generic 
superconductivity, i.e. of a superfluid density ns(T) and a vanishing linear resis- 
tance requires a rigid coupling of the phases of the local droplets. For high- Tc 
materials, the effect of the low carrier concentration on the generic transition 
line has been emphasized recently As illustrated by Fig. la, an estimate of 
the phase fluctuations within the 3D XY-model leads to an upper bound for 
superconductivity 

ksTe = Code- (1-1) 

The rigidity against phase fluctuations is determined by a microscopic cut- 
off length dc |^] for the fluctuations perpendicular to the CM02-planes and by 
the stiffness parameter eg = h"^ /2m* ■ ns{T) = $Q/(47r/ioA^) , which by using 
the flux quantum, $o — /i/2e, can be determined from the penetration depth 
A . In this picture, between Tq and Tc a broad paracoherent regime arises for 
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Figure 1 

Sketches for phase di- 
agrams of high- Tc su- 
perconductors: transi- 
tion temperatures as 
functions of a) hole 
doping H] in zero mag- 
netic field and b) in 
magnetic field near op- 
timal doping. Dashed 
and dotted lines are 
described in the text. 



underdoped materials, 5 < 5 opt , which is related to the widely evidenced pseu- 
dogap phenomena. In real situations, several constraints may enhance the phase 
fluctuations and thus further reduce Tq. At the high temperatures of interest, 
thermal fluctuations decrease the stiffness while finite magnetic flelds - the 
earths field is sufficient for thin films due to their strong demagnetization - in- 
duce a vortex system, which has to be pinned in order to avoid lossy flux creep 
even under applied low current density. The issue of nucleation of generic su- 
perconductivity in the presence of current densities and of vortices pinned by 
various kinds of disorder has been reviewed from the theoretical side |3[ ^, |^. 
Today there is ample experimental evidence |^ for the thermodynamic melt- 
ing line of the Abrikosov vortex lattice Bm(T) in clean materials with different 
anisotropy, 7 = £,ab/£,c (Fig- lb). The location of B,n{T) is described in terms of 
the Lindemann criterium 
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(1.2) 



where aq = ^^q/B is t he mean vortex spacing for _B _L to Cu02 planes. The 
comparison with Eq.(l.l) indicates that the melting is determined by the elastic 



cutoff cj^ao/j given by the small Lindemann number cl- Recent work reveals 
that random point defects give rise to a quasi-lattice (Bragg-glass) at low fields 
{B < Be). Above Bf, entangled vortices (vortexglass, VG) are expected due to the 
increasing repulsive vortex interaction, which drives the vortices into the local 
energy minima set by the fluctuations of the point defects, however, generic 
superconductivity, as conjectured previously |^ for the VG, could not yet be 
distinguished from vortex fluctuations just frozen by uncorrelated disorder 

For correlated disorder, on the other hand, like columnar defects [||piercing 
the entire sample along B, the existence of a generic phase transition Bg(T) 
(Fig. lb) has been evidenced by means of the boson mapping of the vortex lines||^. 
Accordingly, the socalled Bose-glass (BG) phase arises from a localization of 
interacting 2D-bosons in a random potential. It has been argued, that not only 
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defects introduced by heavy ion irradiation but also intrinsic linear structures 
along B arising e.g. by edges, twin colonies or grains of different orientation, or 
screw dislocations may lead to BG-ordcring j^. 

As for the melting line Bm, numerous evidence for generic superconductivity 
mostly for YBCO has been obtained using various experimental techniques. How- 
ever, a detailed understanding of the location of Bg(T) in terms of the pinning 
mechanism in a concrete material, even for well-defined columnar defects, has 
not yet been achieved. In this work, we investigate Bg(T) in c-axis oriented, gran- 
ular YBCO films with large critical current density, jc{T ^ 0) = 2 • lO^A/cm^ 
grown for SQUID fabrication In addition to the intrinsic disorder, we study 
the influence of columnar defects of various density in order to compare intrinsic 
with correlated pinning. We also investigate the singularities of the fluctuation 
conductivity near Tg[B) in order to shed more light into the nucleation process 
of generic superconductivity in real high- Tc materials. 



2 Inductive Probe of Generic Superconductivity 

Along with direct, four terminal I/V - and ac-probes of the superconducting 
transition, for the brittle cuprates inductive techniques became very popular, 
because they avoid contacts and turn out to be high sensitive to the high con- 
ductances near of fllms or bulk samples. Most of them employ irreversibility 
effects in the magnetization, the screening, the onsets of nonlinearity or of higher 
order harmonics in the ac-susceptibility x': ^^nd maxima in the absorption x" to 
detect Tg{B). Similar as with the voltage criterion in I/V-curves, commonly used 
to define the transition temperature, also these methods are based on at least 
one extrinsic condition, like frequency, measuring time, excitation amplitude, 
sample size, or simply the detection limit. 

Here we use the frequency dependence of the magnetic susceptibility measured 
between 30mHz and 30MHz with amplitudes bo{uj) from lO""* to 10~^mT to 
keep within linear response and to optimize the sensitivity. Fig. 2a illustrates 
that the current density and hence the conductivity cr(a;, T) in the Cu02 planes 
is probed. For long cylinders, where demagnetizations are negligible, xi^) can 
directly be obtained by integrating the solutions of Helmholtz' diffusion equation, 

h = {plq(t)^^ l\h as used e.g. in Ref. The self-field effects produced by thin 
films, lead to a non-local diffusion equation, which has been solved numerically 
for all practical sample geometries and field configurations by E.H. Brandt p^ . 
The results are presented in the form 
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Figure 2 a) Contactless technique for measuring liigii resistivity parallel to the 
Cu02 planes; b) dc-resistivity measured in zero field near Tc on a heavily twinned 
Y Ba2Cu^07-s crystal (Lz = 800^m) and a thin film (Lz = 0.2/im). Indicated 
are power law singularities, p ~ {T /Tc — ly , with s=7.7 fitting the data and 
s — 1/3 (dotted) predicted for clean superconductors. 



where the coefficients c„ and A„ become independent on sample size for « 
R. The complex conductivity in k'^{uj) = iloR"^ fio{a' (uj) + ia"{ui)) is extracted 
from x{^) using a inversion routine. For the present precision, Sx/x ~ 10""*, n = 
30 terms are sufficient, which for our case of circular discs have been tabulated 
in Ref. jl^, where also more details are given. 

As example, we compare in Fig. 2b the zero-field dc-resistivities of the present 
thin film, = 0.2/im, with that of a heavily twinned crystal, = 800/im, 
p| . Despite the different T^s, 88K and 91K, the results fall on the same power 
law, which is characteristic for a continuous phase transition. This also indicates 
some universality, which, however, deviates significantly from the prediction of 
models for dean materials: the mean-field, s = 1/2, and the 3D-XY with s = 
1/3. 

The ac-conductivity near Tc, as shown for a' in Fig. 3, reveals no clear signature 
of a sharp transition. The values of a are rather high as compared to the Drude 
background l]!^ . They arise from the low frequencies which indicates strong 
superconducting fluctuations, 

a(c., T) = -i- I ^ + r„ 1 + a/, (^, T). (2.4) 



In the first, the 2-fluid term, — m* / fiQ^'^i^e denotes the plasma wavelength 
and = 2kTc/fi = 2.5 ■ lO^^s"^ the large scattering rate of the quasiparti- 
cles. In the absence of fluctuations, the phase angle a" /a' provides the strongest 
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Figure 3 Real part of the 
dynamic conductivity multi- 
plied by frequency near the 
transition Tc = 90.957^ 
(s. Fig. 4) to generic (sta- 
ble) superconductivity of a 
YBa2Cu307-s film in zero 
magnetic field. 



signal at Tc by jumping from zero to Us/uJTn. Due to Ufi this jump is smeared 
by screening and lossy contributions appearing above and below Tc^ respectively. 
They produce a positive Drude-like frequency gradient d{a" / a') / duo > above 
Tc and the inverse behaviour, d{(j" /a')/du) < below Tc, so that for not too 
high frequencies, wr^ << 1, at Tc a crossing of all a" /a' curves can be expected. 
In fact, plotting a" /a' in Fig. 4a, a well-defined crossing temperature is realized. 
This signature for nucleation of generic superconductivity in both, the Meissner 
and vortex states as well was pointed out by D. Fisher et al. p] and is uti- 




Figure 4 Temperature dependence of the phase angle of the dynamic conduc- 
tivity, a' — ia" a) in zero magnetic field determined from the data of Fig. 3 and 
b) at i3 = 19r for the same film. Insets display the modulus |cr|. 
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lized here to trace the transition hne. As demonstrated by Fig. 4b, the crossing 
point persists up to the maximum available field of 19T. It is, of course, shifted 
downward in temperature but the height at T = T^, a" /a' — 3.7(2), remained 
unchanged. Both features will be discussed in section 4. 



3 Fluctuation Dynamics and 'Generic' Transition Lines 

Since the spatial extension of the thermal fluctuations diverges at the continu- 
ous transition, ^(T) = ^(0)|1 — T/Tgl"'', all macroscopic quantities at or near to 
equilibrium should be invariant against scaling by some power of this correlation 
length ^, like, for example, the relaxation time of the long-wavelength fluctua- 
tions, T ~ Scaling has been proposed for the phase o'"/cr' |^ by using the 
scaled frequency, lot. They collapse onto two branches of the data for above and 
below Tg is demonstrated in Fig. 5 for cr(w) measured in zero field and in B=1T. 
In addition to the scaling property we note that the presence of the high vortex 
density did not change the shape of the two scaling functions, a" /a' = P±{lot). 
Moreover, the critical slowing down of the relaxation rate 

T-\T,B)^T^H'^^ (3.5) 




Figure 5 Dynamical scaling of the phase of the dynamic conductivity for the 
pristine film a) at zero field, where Tg{0) = Tc = 90.85-?^ and b) at B = ITesla, 
where Tg = 86AK. 
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can be described by a single empirical form, if the reduced temperature difference 
t{T, B) = \T- Tg{B)\/{To - T) is introduced suggested in refs. 0. By this 
choice, both the exponent vl — 9.5(3) and also the amplitude t^^ = 0.8- lO^^s^^ 
remain the same on either side of Tg and for all applied fields with Tq = Tg{Q) + 
O.SOif as the 'mean-field' temperature. 

All these features provide strong evidence that the onset of generic supercon- 
ductivity along Tg{B) is triggered by a common mechanism independent of the 
field induced vortex density. Since the phase angle at Tg is given by 

LPg ^ arctan ^ ^ ^(1 - (3.6) 

this also explains the field independence of ipg, i.e. z=5.7(3) in Fig. 4. Moreover, 
the Kronig-Kramers relations imply a power law at T — Tg, a ^ (iwtq)^/^"^ 
and a scaling behaviour of the modulus, if \a{uj, T)\ is scaled by the appropriate 
static limits, which are here ao{T > Tg{B)) ~ t~'^ with s=7.7 (see Fig. 2b) and 
(To(T' < Tg) ~ /u) with V = 1.7(f). The unversality of the critical exponents for 
the vortex state has been emphasized in Ref. ||l^ based on numerous reports for 
thin films, which now include the zero-field limit for YBCO and, amazingly, 
also ac-data from 150nm thin indium film Jl6| . A change of the universality 
class, however, was realized after implementation of columnar pinning centers 
by heavy-ion irradiation with a dose nc = 2 • lO^^cm^^ [0. Here we extend this 
study to higher and lower matching fields, Brj, — riccj^o. 

At first, let us look at the critical exponents v and z, evaluated from ipg 



(Eq.(p!q)) and from the slowing down of (Eq.(3.5)). The results shown in 
Fig. 6 clearly demonstrate significant changes through the defects, except for 
B < lOmT to be considered elsewhere IfSl. We find that the irradiation reduces 




Figure 6 Field depen- 
dence of the critical expo- 
nents 1/ and z determining 
the divergences of the corre- 
lation length and the relax- 
ation times associated with 
the transition to generic su- 
perconductivity. 
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Figure 7 Effects of mag- 
netic field B\\c and colum- 
nar defects B^\\c on the 
transition temperatures Tg 
to generic superconductiv- 
ity. Inset: influence of on 
Tg in zero magnetic field. 
Full lines represent power 
laws being discussed in the 
text. 



z, but since the exponent of , j^z=9.5(3), remains unchanged, v rises by the 
same amount as z decreases. Upon filling the defects with vortices, the exponents 
change rather sharply at some field < B^, where they become identical to 
those of the pristine film. We qualify B^ as an effective matching field being 
somewhat smaller than B^. Similar differences between B^ and the actual defect 
densities have been reported by direct images of the pins . 

Another dramatic effect of the irradiation occurs in the generic transition 
temperatures shown in Fig. 7. We note four significant features : (i) the transition 
lines of the pristine film and the film with small pin density, B^ = IT, fall on 
the same power law, indicated by the thick line: 

Bg{T,B^) = Bg{B^){Tg{B = 0,B^)/T-lf (3.7) 

with (3 = 1.33(3) and Bg = 50(2)T. (ii) At larger doses, the transition temper- 
atures are drastically reduced in contrast to what one would naively expect for 
strong correlated pinning. This can be well described by a fit to Eq.(p77|) with a 
greater exponent /3=4.0(3), indicated by thin lines in Fig. 7. (iii) If there are more 
vortex lines than pinning rods, B > B^, the transition lines approaches that of 
the pristine film, which indicates that the pinning by the rods becomes ineffec- 
tive, (iv) As shown by the inset to Fig. 7, the zero-field transition Tg{B = 0, i?^) 
is suppressed by irradiation, which is a well-known feature but has not yet been 



explained. Accordingly,this reduction of Tc obeys the power law, Eq.(3.7), with 
/3—1.33{7) and the same amplitude i3g=50T. This implies that the columnar 
defects exactly act like vortices in reducing the transition temperature. With 
slightly larger amplitude the same law holds for the transitions in irradiated 



YBCO-crystals extracted from Ref. [20 
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4 Phase Fluctuations in the Presence of Disorder 



We start the discussion with two interesting observations made on the pris- 
tine film. The first is that the transition line Bg{T), Eq.(3.7), turns out to follow 
rather closely the melting line predicted by Eq. (^]^) : inserting the measured pen- 
etration depth, for which in the present range of temperatures we take A((r) = 
140nm/ (l-r/re)i/3 wegetforthe stiffness eo(T) = 980ir(l-T/Tc)2/3. 

Then Eq.( |1.2D yields /3 = 4/3, in perfect accord with the fitted temperature vari- 
ation of the data in Fig. 6 in fields up 20T. Moreover, by comparing the predicted 
amplitude Bg — 0o[c|eo(O)/7rc]^ with the measured value of 50(2)T and using 
7 = 5 for YBCO we obtain for the Lindemann number = 0.26(1). This value 
agrees surprisingly well with very recent numerical simulations: (i) within the 
Bose-model for the vortex lines, cl = 0.25 and (ii) also with cl — 0.24 
from the anisotropic 3D-XY model for a London superconductor, A/^ >> 1 p^ . 
However, unlike the Bose model, the 3D-XY simulations also comprise the zero 
field limit where no external flux lines present. They indicate that the transitions 
in both zero and finite field are driven by the Onsager unbinding of thermally 
excited, closed vortex loops. The result ~ (1 — T/Tc)"-^'' [|2^ is fully consis- 
tent with the measured penetration depth by which we explained the exponent 
of the generic line (3 = 4/3. For B = 0, Williams |^ arrived at essentially the 
same variation of tt-s and interpreted it by a blowout at Tc = Tg{0) due to a 
percolation of vortex loops determined by the sample size R. At finite B, the 
externally induced vortex-lines appear to mediate a percolation already between 
small loops, thus reducing the transition temperature Tm(B). 

With this in mind it is suggestive also to attribute the reduction of the zero- 
field transition in the co/Mmnar defected films (see inset to Fig. 7) to the thermally 
excited loops. We assume, that the rods act as nucleation centers for the loops 
which save there some nucleation energy. Then the same amplitude of 50T de- 
scribing the vortex-induced and rod-induced reduction of Tc indicates that the 
saving is about the same. The suppression oiTc{B^) is described by the dephas- 
ing criterion, Eq.(1.2), if in the phase cutoff, dc = cj^Ur/j, the mean distance 
between the rods Oj. instead of the mean vortex distance oq is used. 

In finite applied fields, the transition temperature of the columnar defected 
films would stick at Tc{B^) up to _B = if all induced vortex lines would localize 
at the rods. Obviously this is not the case so that we consider the possibility of 
thermal depinning. Here we base our discussion on the Bose model since it was 
argued recently |Q , that for S > some results like cl turn out to be the same in 
the vortex-loop and the Bose-model. Accordingly, we use for the depinning from 
rods a phase cutoff at dr = tq /tt^ |9| . By inserting this and the mean pin radius 
in our films, tq = 3.5nm, into Eq . (|l . l| ) , one finds Tdp = O.STc = 72/^, so that in 
fact the delocalization is occuring in the temperature range of interest, see Fig. 7. 
However, since the BG-work treats the thermally smeared pin potential only as 
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a perturbation of the vortex- vortex interaction for B > Bd , the transition line 
is always shifted to above the melting hne (see Fig. lb), while our data approach 
Bg{T, Btp) ~ B,n{T) from below. We believe that our results should be analyzed 
in a more consistent picture, which takes into account the phase breaking by 
thermal loops in the presence of columnar pins and externally induced vortices 
. One step into this direction has been made by Wallin and Girvin |p5| who 
obtained from simulations based on isotropic vortex loops v = 1.0(1) and z = 
6.0(5) . Interestingly, these values are fairly close to our low-field results, but 
they clearly fail above lOmT, i.e. in the regime where Bg{T) rises steeply with 
/3 = 4 in Eq.(3.7). On the other hand, the agreement with the exponents of 
the pristine film indicates that there intrinsic linear defects nucleate the generic 
transition. 

The most obvious consequence of the increase of the exponents over z = 3/2 
p3[ and v — 2/3 |22j for a clean 3D-XY system is the broadening of the fluc- 
tuation regime. This is illustrated for the dc resistance by Fig. 2b for which one 
expects po ~ (r — TcY'^^^^K Qualitatively spoken this broadening may arise 
from an increased nucleation of vortex loops below Tg and from a reduction of 
these phase fluctuations above Tg by the intrinsic disorder. Similar effects on the 
correlation length i.e. on v have been observed earlier on structures consisting 
of weakly Josephson-coupled Nb-grains To best of our knowledge, the in- 
terplay between the thermally excited loops, field induced vortices, and the real 
disorder in high- Tc materials has not yet been worked out. 

To conclude, we provided evidence that thermally excited vortex loops de- 
termine the nucleation of generic superconductivity of YBCO-films in zero and 
finite magnetic fields. This follows from the identical scaling behaviour of the dy- 
namic conductivity for B = and B > and the quantitative description of the 
generic transition line by the 3D-XY, i.e. vortex-loop model. Additional support 
for the dominance of these phase fluctuations comes from the observation ||l^ 
that the relevant temperature variable, which described the critical behaviour 
up to vortex densities B = 19 T, is given by t = |T - Tg{B)\/{To - T) (Eq.(U). 
Such behaviour emerges from the temperature variation of the coupling energy 
between the phases in the 3D-XY Hamiltonian ||2^, |2^ . 
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